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We study spreading on networks where the contact dynamics between the nodes is governed by 
a random process and where the inter-contact time distribution may differ from the exponential. 
We consider a process of imperfect spreading, where transmission is successful with a determined 
probability at each contact. We first derive an expression for the inter-success time distribution, 
determining the speed of the propagation, and then focus on a problem related to epidemic spreading, 
by estimating the epidemic threshold in a system where nodes remain infectious during a finite, 
random period of time. Finally, we discuss the implications of our work to design an efficient 
strategy to enhance spreading on temporal networks. 


I. INTRODUCTION 

When modelling spreading in systems, it is often as¬ 
sumed that an entity, e.g. a virus or information, per¬ 
forms jumps on a static underlying structure represented 
by a network [1]. In its simple form, the network is as¬ 
sumed to be undirected and unweighted, and it can be 
made more elaborate by incorporating additional infor¬ 
mation about the direction of the edges and their weight. 
In this framework, spreading is described by linear dif¬ 
ferential equations, where the density on nodes evolves 
according to the flux induced by neighbouring nodes, and 
whose solutions is determined by the spectral properties 
of a matrix encoding the structure of the network, of¬ 
ten the Laplacian matrix. In a large number of systems, 
however, this modelling approach is not justified because 
events taking place on the nodes exhibit complex tem¬ 
poral patterns, and the underlying structure has to be 
considered as a temporal network [2]. For instance, in 
systems as diverse as mobile phone communication, email 
checking and brain activity, nodes and links are not per¬ 
manently active but are punctually activated over time, 
and their dynamics tends to deviate strongly from that 
of a Poisson process. 

An important research question associated to the tem¬ 
porality of the network is to understand its impact on 
the speed and reach of spreading IMS]. For instance, 
the presence of correlations between activations of neigh¬ 
bouring edges may either slow-down or accelerate spread¬ 
ing by inducing non-random pathways, depending on the 
type of correlations mm- Another important mechanism 
is associated to the burstiness of the temporal processes, 
as two consecutive activations of a link or a node tend 
to present a broad distribution of inter-contact times [8] . 
The latter mechanism is the subject of this paper. Our 
main contribution is to consider a model where spread¬ 
ing is not always successful on an edge when it appears, 
so that it takes place only after a random number of tri¬ 
als. This process allows to introduce, and to tune, differ¬ 
ent time scales in the system, one associated to network 
evolution, through its inter-contact times, and the other 
to spreading. As we show analytically and numerically, 
imperfect spreading leads to a generalisation of the so- 
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FIG. 1. Illustration of the sequence of activation times of an 
edge. The inter-contact times ti are randomly taken from the 
distribution /(r). The infection of one end of the edge takes 
place at a random time, illustrated by a star. 

called bus paradox, and has interesting implications on 
the properties of epidemic spreading on networks and, in 
particular, on how to optimise the spreading of a spread¬ 
ing process. 


II. SPREADING ON TEMPORAL NETWORKS 
A. Estimating the inter-success time 

The spreading model is defined as follows. Let us con¬ 
sider a temporal network where the activation of edges 
is governed by a renewal process, with inter-contact dis¬ 
tribution /(r). It is further assumed that edges have a 
vanishing duration, and that the random processes as¬ 
sociated to the different edges are independent, spread¬ 
ing takes place from node to node. The diffusing entity, 
called the virus from now on for the sake of simplicity, 
sits on a node until an edge appears. The virus then 
invades the neighbouring node with a probability p. In 
order to describe the random process, it is crucial to es¬ 
timate the waiting time before a certain edge is invaded, 
which we call the inter-success time. So far, this model is 
essentially equivalent to an SI model, which we assume 
to take place on a tree in order to avoid non-linear ef¬ 
fects. In situations when p = 1, the relation between the 
inter-contact distribution and the inter-success distribu¬ 
tion is well known to lead to the so-called bus or 
inspection paradox, as the mean inter-success time (t) is 
given by (r^) /2 (r), where (r) and are respectively 
the first and second moments of the inter-contact time 
distribution. In situations when the inter-contact distri- 
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bution presents a broad tail, the average success time can 
therefore become arbitrarily large as compared to the av¬ 
erage contact time. This result, which seems paradoxical 
at first glance, originates from the fact that the proba¬ 
bility for the virus to arrive during a inter-contact time 
interval is proportional to its duration. Because the acti¬ 
vations of different edges are independent, it can indeed 
be assumed that a virus arrives at a random time within 
an inter-contact interval (see Fig. 1). Based on these 
arguments, one finds that the probability of a success in 
a first attempt at time t is 

= /Wdr. (1) 

Let us now turn to the case of an arbitrary value of p. 
From the previous expression, the inter-success probabil¬ 
ity after an arbitrary number of attempts is given 

by 

oo 

V’W = ^p(l-p)''P(fc + l,i), (2) 

k=0 


where p{l—p)^ is the probability of a success after /c + 1 
trials, 

p{k+i,t) = iQ*r^)it), ( 3 ) 

where * stands for the convolution product and P{k,t) is 
the inter-success probability in k trials. This expression 
takes a simple form in the associated Laplace space 

Pik + l,t) = Qis)p{s), (4) 


where the upper tilde stands for the Laplace transform. 
Using the properties of the small s expansion of each 
distribution, 

^{s) = l-s{t) + ^ {t^) + 0{s^) 

/(s) = l-s(T) + y(T2) + £)(/), (5) 

one directly find a relations between their moments with, 
in particular for the first moment. 
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and we recover the expressions of the standard bus para¬ 
dox when setting p = 1. 

Decreasing p systematically increases the average inter¬ 
success time, as expected because more and more trials 
are required for the spreading to actually take place. In 
order to account for this trivial effect and to properly 
identify how the shape of the inter-contact distribution 
affects (t) as a function of p, we focus on the normalised 
average success time f, a standard measure for the bursti- 
ness of a process HHg, defined as 


T = 


it) 

{t) /p' 


( 7 ) 



FIG. 2. Proportion of infected population (dashed green) and 
reproduction number (solid blue) for a tree-like network as a 
function of p. (n) P = 1 indicates the epidemic threshold. 
The parameters of the model are (n) = 2, r(t) — — and 

the inter-contact distribution is a Gamma distribution, with 
parameters a = 0.5 and 5=1. 


where (r) /p is the average success time in the case of a 
Poisson process. In the latter case, the inter-contact time 
is exponential and = 2 (r)^. One directly finds 

an expression taking the value 1 in the case of a Poisson 
processes, as expected, and depending linearly on the so- 
called burstiness coefficient 


P = 


2{rf 


- 1 , 


which takes a positive value in fat-tailed distributions, 
as the ones observed in a majority of empirical systems. 
This result clearly shows that burstiness tends to slow 
down dynamics, but that its impact becomes less and 
less important as the probability of success p decreases 
and a higher number of attempts is necessary for the 
virus to spread. 


B. Epidemic threshold and spreading efficiency 

The dynamics described so far allow to estimate the 
spread of infection in situations when nodes do not re¬ 
cover from the disease. In a majority of practical situa¬ 
tions, however, nodes remain infected only during a finite 
time before recovering. In that case, the propensity of 
the virus to invade the system is mainly governed by the 
connectivity of the network, e.g. number of contacts per 
infected user, and its transmissibility P, also called infec- 
tivity, defined as the probability that the virus spreads to 
an available neighbour before the infected node recovers 
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FIG. 3. Isolines for the transmissibility P in the {{r),p) plane, 
for Gamma distributions with different shape parameters a 
and r{t) = — 1). In the case of a Poisson process, isolines 

are linear, indicating a dependency of P on {r)/p. Dashed 
lines indicate the linear relationship of a Poisson process. De¬ 
viations from a = 1 alter the convexity of the isolines. 


and the virus becomes de-activated. Transmissibility is 
defined as 

P+OO r+oo 

P= / '0(t) / r(r)drdt, (9) 

Jo Jt 

where, as before, '0(t) is the probability of a successful 
infection at time t, and r(r) is the probability that the 
infected node recovers at an ulterior time r > t. This ex¬ 
pression clearly shows the importance of the competition 
between two temporal processes. In the case of tree-like 
networks, where all nodes have the same transmissibility 
P, it is straightforward to show that the basic reproduc¬ 
tion number, defined as the average number of additional 
people that a person infects before recovering, is given by 

Ro=nn), ( 10 ) 

where (n) is the expected number of susceptible neigh¬ 
bours of an infected node. The epidemic threshold is de¬ 
fined by the condition Rq = 1 separating between grow¬ 
ing and decreasing spreading. The epidemic threshold is 
thus reduced either by reducing the transmissibility or 
(n). This result is valid only when the network has a 
tree-like structure, which is valid for a majority of ran¬ 
dom network models below the epidemic threshold m- 
In order to illustrate these results and to perform 
numerical tests in the following, we consider an inter¬ 
contact distribution given by a Gamma distribution 


/(t;a,6) = -—^ 

^ ^ h^r{a) 

with scale parameter b and shape parameter a. The fam¬ 
ily of distributions has the advantage of including the 
exponential distribution, for a = 1, and to produce dis¬ 
tributions with a tuneable variance by changing a. One 
finds 



FIG. 4. Transmissibility P vs probability p of success for fixed 
values for fixed values of {r)/p, with r{t) = S{t — 1). In the 
case of Poisson processes, with a = 1, the system does not 
exhibit a dependency on p. When the system is bursty, in 
contrast (a < 1), increasing p decreases the transmissibility 
of the process. 


(r) = ab 

(11) 

a 

(12) 


Fixing (r), one thus finds that the variance increases by 
decreasing a. 

Let us also note that transmissibility takes a particu¬ 
larly simple expression when / is an exponential distri¬ 
bution 1 /so that 

P = 1 - (13) 


when the recovery time is a constant tc, r{t) = 5{t — tc) 
and 


ptc/ {t) 

ptc! {t) +1 


(14) 


when the recovery distribution is an exponential with 
rate l/tc- As a first check, one shows the accuracy of 
Eq. to determine the epidemic threshold on a tree by 
numerical simulations in Fig. 

Let us now consider the impact of the shape of the 
Gamma distribution, calibrated by a, on transmissibil¬ 
ity, by using the exponential case a = 1 as a baseline. 
Numerical solutions show in Fig. the dependency of 
transmissibility on (r) and p, for three values of a (the 
value of b is thus determined because of Eq. 11). Simu¬ 
lations are all performed with r{t) = 6{t — ijbut simi¬ 
lar results are obtained for exponential distributions for 
the recovery time. One observes that transmissibility de¬ 
pends on the ratio {r)/p when a = 1, but that the sys¬ 
tem exhibits deviations to this linear relationship when 
the dynamics deviates from a Poisson process. In order 
to quantify this deviation, we explore the dependency of 
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FIG. 5. Transmissibility as a function of the recovery time tc, 
for r{t) = S{t — tc) and fixed values of {t)/p = 1. The dynam¬ 
ics is more and more affected by the shape of the distribution 
for larger values of p. 

P on the success probability p, for fixed values of the pa¬ 
rameter {r)/p. The latter is simply a naive estimation of 
the time of success, obtained by multiplying the average 
time between two contacts (r) and the expected number 
of contacts before a success takes place 1/p. Numeri¬ 
cal results clearly show, in Fig. that transmissibility 
becomes less and less efficient, as compared to the Pois¬ 
son case, for increasing values of p, in agreement with 
the observation that the inter-success time is more and 
more affected by the bus paradox for increasing values of 
p. The mechanism behind this effect is also apparent in 
Fig. where, again for a fixed value of {r)/p = 1, trans¬ 
missibility is clearly less efficient in the case of bursty 
dynamics, for larger values of p. 

This result has important practical implications, for in¬ 
stance in online marketing, where user activity is known 
to be bursty m- Indeed, let us consider a marketing 
agency having to decide on a strategy in order to opti¬ 
mise its viral impact. The impact of its campaign, eval¬ 
uated by its transmissibility, increases only sub-linearly 
with its quality, evaluated by its success probability p. 
This observation therefore needs to be considered when 
devising a strategy to balance the gains and the cost of 
increasing the viral potential of an ad. 

III. CONCLUSION 

The main purpose of this paper was to study spread¬ 
ing on temporal networks, where the temporality of edge 


activations is modelled as a stochastic process and where 
the inter-contact time distribution takes an arbitrary ex¬ 
pression. In practice, we have mainly been interested 
in situations when the inter-contact time distribution 
presents a fat tail, as observed in a variety of social net¬ 
works. In contrast with a majority of previous works, we 
have considered a model where the transmission of the 
diffusive entity is only successful with a certain probabil¬ 
ity at each contact. Our main results are twofold. First, 
we have derived an analytical expression for the average 
inter-success time. We have shown that burstiness plays 
a more and more important role, associated to a slowing 
down of the dynamics by the so-called bus paradox, when 
the probability of success is increased. In the limit of a 
small probability of success, the shape of the distribution 
ceases to play a role, and only its average determines the 
speed of spreading. Second, we have turned to numerical 
computations in order to calculate the epidemic thresh¬ 
old for a dynamics where nodes remain infectious during 
a random period of time before recovering. Our results 
confirm that the bus paradox hinders the spreading of the 
process when the probability of success is increased. As 
we discussed, the results have implications for the design 
of efficient marketing campaign. In particular, it suggests 
that an efficient strategy should aim at properly predict¬ 
ing the future activation times of a target user, in order 
to minimise the time between his infection and his first 
contact, and therefore the impact of the bus paradox. In 
this work, we have considered locally tree-like networks, 
as often assumed when studying spreading processes. An 
interesting line of future research would be to incorporate 
more complex, arbitrary network structures, where cycles 
and communities play a role US]. 
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